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I suggest that electron localization by random flux and passive transport in quenched velocity fields 
in two dimensions be studied as perturbations of the simple operator IC — A • V, with incompressible 
velocity field/vector potential A = V x = {—dy, dx)4'- This operator has an infinitely degenerate 
subspace of zero energy eigenstates, arising from incompressibility, that are extended ior generic 0(x) 
and are expected to remain so under perturbation. I propose that an anomaly accounts qualitatively 
for properties of the spectrum and eigenstates of IC and its perturbations. 
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In two dimensions, all eigenstates of the hermitian op- 
erator Hv = — V^ + V{x) are localized for random V, 
no matter how small the randomness (Anderson localiza- 
tion). Over the years it has emerged that certain kinds 
of fields can suppress localization. In particular, a strong 
uniform magnetic field B = VxA = dxAy—dyAx produces 
a dramatic change in the character of the spectrum of 
TCb = (— iV — A)^/2m, yielding Landau levels consisting 
of degenerate extended eigenstates that remain extended 
on perturbation by random V. 

The underlying physics of this phenomenon, the quan- 
tum Hall effect, can be thought of as arising from in- 
compressibility. In brief, taking ?7i ^ in the La- 
grangian Cb = rni^/2 + A-x, one obtains the action 
fdtA-Ji. = J^cPaB. This action is manifestly invariant 
under area-preserving deformations of il. The invariance 
persists after quantization - albeit in a different form - 
and is the origin of 'quantum incompressibility' O. 

This version of incompressibilty is expected to apply 
when there are spatial fluctuations of magnetic field, pro- 
vided the mean magnetic field is non-vanishing. When 
the magnetic field has zero mean one enters a differ- 
ent regime that is perturbatively inaccessible from Lan- 
dau levels; rather, a different kind of incompressibility 
may provide a more suitable starting point for pertur- 
bative analysis. Such a regime can be achieved in two- 
dimensional systems with quenched spatial disorder, such 
as the random flux (RF) model Hrf = (-iV + A/2)^ 
(in gauge V • A=0 for our purposes), and the non- 
hermitian passive scalar advection-diffusion model (PS) 
'Hps = V^ + V • (v-), with velocity field v = A + u and 
vorticity w = V x A, where V- A = and V x u = 0. 

Common to both these models is a component of the 
form IC = A-V; for example, 7^^ = VV2m + A- V ^ /C 
as TO — > CX3. As shown in this Letter, incompressibility 
determines the character of eigenstates of K. even when 
B (or w) = — V^(/) is not uniform. Models containing K, 
are ubiquitous, and basic features of the spectrum and 
eigenstates of K, can be expected to survive perturbation 
by — V^ and A^ or V provided the perturbations are 
sufficiently weak. 

The quenched random fields of the models fall into two 



classes. The vorticity w = Vx A in Tips is, like magnetic 
field, a pseudoscalar quantity, in contrast to the scalar V- 
u. Similarly, the randomness in Tirf can be decomposed 
into parity-odd and -even contributions iA- V and A^/4. 
This distinction governs the contribution to scattering, 
and suggests the study of models, which I call RFI and 
PSI, that retain only the parity-odd random components 
of RF and PS respectively, with equations of motion: 

(^)9tV = V•J^-V■{(-)VV^ + (^)A^} (1) 

(The factors of i and — 1 should be chosen appropriately 
for each model, but with respect to the symmetry ad- 
dressed in this paper these factors are inconsequential .) 

The divergence form of (1) shows explicitly that 
RFI and PSI conserve three currents |^: the ordi- 
nary probability-current, conserved also when parity- 
even fields are present, but in addition the ■(/'-current J^ 
generated by ip ^ ip + C and its transpose f/^^-current. 
When second quantized, the currents correspond to three 
non-commuting operators: the number operator Nq = 
^^qci, and the V^-current operators cj = X^j*^! ^^'^ 
Co = '^iCi, the sums being over positions in space. All 
three operators commute with V^ and parity-odd opera- 
tors, but only A'o commutes with potential disorder. The 
other two fiuxes embody a dynamical constraint B on re- 
laxation of both the real and imaginary parts of ip: when 
it is satisfied, Keip- and Im-f/'-fluctuations leaving one re- 
gion must be exactly compensated somewhere else; V,u, 
and A^/4 violate the constraint by introducing uncom- 
pensated sources and sinks of Re'0 and Im-0. When the 
current J^ and its transpose are conserved, I say that 
the coupling to the random field is incompressible M . 

In obtaining RFI from RF by neglecting A^ , I have fol- 
lowed ref . [6] , where standard methods are then used to 
map RFI to a unitary non- linear sigma model (NLSM), 
for which a/Z eigenfunctions are localized. From the local- 
ization of all eigenfunctions of RFI, it is deduced in ref. [6] 
that all eigenfunctions of RF are localized. In contrast, I 
show below that an infinite number of eigenstates of RFI 
are extended, which suggests that when A^ can be prop- 
erly treated as a perturbation to RFI, the corresponding 



eigenfunctions of RF will be extended. The conservation 
laws (1) should be reflected in the correct NLSM for in- 
compressible disorder, but have no counterparts in the 
unitary NLSM, indicating that the conclusions of rcf. [6] 
may not be universally valid. 

Generically, iA-V does not commute with — V^, and 
the spectra of these operators differ qualitatively. If 
HC^p = Eip, then i/CV'^ = —Eij}'^ and i/CV" = nEip^ , so 
the spectrum is symmetric and unbounded. It follows 
that at zero temperature, whereas the k = mode of 
a free particle is buried at the center of the fermi disc, 
for jA • V this mode sits instead at the center of the 
spectrum, yielding a dirac sea rather than a fermi sea. 

To further study RFI, I at first examine the spectrum 
and eigenfunctions of /C alone, and only later restore 
— V'^. Eigenfunctions oi ilC = iV x </) • V lie on char- 
acteristics: the trajectories of x = A, or equivalently, 
the streamlines (equipotentials) of (f). The symmetry of 
/C is elucidated by observing that ICip is just a Poisson 
bracket, if (j) is viewed as a Hamiltonian and y as momen- 
tum conjugate to x: ICip = {ip, (j)}pB =dxipdy4' ~ dy4'dx4>^ 
so that the time evolution of ij) is given by = di/j/dt = 
dti^ + {(j),'4'}pB- Incompressible reparameterizations are 
the canonical transformations of (a;,y). /C is covariant 
under this non-abelian group of transformations g . 

Under a subgroup of the canonical transformations, 
/C is invariant. This subgroup represents translations 
along the streamlines of 0; in infinitesimal form (77 <C 1), 
X = x' -I-77V X f{4>), so that (/)(x) = 0(x'). These transla- 
tions are generated by f{4>), which satisfies Vx(j)-Vf{4>) = 
and represents an infinite set of conserved charges, 
one for each streamline. The Lagrangian correspond- 
ing to (j){x,p) is not invariant under these translations, 
but varies under canonical transformation by a total time 
derivative |8|, which can lead to an anomaly |H|. 

Time-independent canonical transforma- 

tions act isospectrally on /C; it is instructive to rewrite 
the Hamiltonian in action-angle variables. Take s to be 
arc-length along (j){x.) — 0o and set x = Jq d.s' /\\7(j)\^ 
where x(s = 0) is an arbitrarily chosen point on the 
streamline. Defining the period t{(J3q) — §, dx, one ob- 
tains as variables the area enclosed by the streamline, 
J{4>q) — Id, ydx (action), and w — x/t (angle). The 

eigenstates of iK. = i(r(0o)) dw on cquipotcntial 00 
are 
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with eigenvalues Em.ipo — '2.'Km/T{4'{)). In this local pa- 
rameterization, the form of the Hamiltonian is (half) a 
chiral fermion, yielding the dirac sea obtained above. 

It is helpful to examine the simpler case when is 
a function of only the coordinate x, where x lies on a 
torus r^ of side L\ (j){x + L) = 4>{x). The spectrum is 
then {{2nm/L){dx(j){x)}xo '■ < xq < L}. Zero modes 
are functions independent of y, and the fourier modes 
exp (z27ma;/L), n integer, constitute a basis of extended 
eigenstates a.t E = Q. Each non-zero eigenstate S{x — 



xa)f{x)exp{i2'Kmy/L), for integer rre, corresponds to a 
one-dimensional wave-function, confined transversely to 
Xf) and extended longitudinally in y around the torus. 

To see that the non-zero eigenstates are extended in 
the y direction, one applies a phase twist a at y = 0. A 
unitary transformation by expiax/L then yields iJCa — 
{dx(l)){idy -t- a/L). As a varies from to 27r, the eigen- 
values for any fixed xq shift up or down an energy level 
depending on the sign of dx4>, whereas the wavefunctions 
remain constant. This sensitivity of eigenvalues to per- 
turbation at the boundary is known as 'spectral flow,' 
and is the hallmark of an anomaly |^ . 

So long as (j) is independent of y, the model can be 
directly rewritten in fourier basis on an iV x A^ lattice. 
Of the 7V^ modes, iV^ — N are extended in the y direction 
and localized in the x direction, and the remaining N zero 
modes, the m — Q fourier mode on each streamline, are 
extended in both the x and y directions. 

For generic (J3{x,y), the non-zero eigenstates remain 
one-dimensional and conflned to the streamlines of (/'(x); 
however, they no longer all traverse the torus. One ex- 
pects generically two sets of streamlines with opposite 
orientation that circulate one of the holes in the torus, 
but most eigenstates close without circulating the torus; 
they are localized. The zero modes are once again com- 
posed of functions that are constant along the stream- 
lines, but have arbitrary variation perpendicular to the 
streamlines; e'^^™'^ for integer n constitute a suitable ba- 
sis of extended eigenstates at E = 0. Eigenstates on 
streamlines that circulate the torus will once again shift 
by an energy level as the boundary phase twist is in- 
creased, but the localized closed loops will undergo spec- 
tral flow only if a phase twist is imposed somewhere along 
the loop. 

The eigenstates may be further characterized by defin- 
ing translation operators T^ — e''^'' for fi — w, J cor- 
responding to the pseudomomentum operators id^. T^ 
translates wave functions parallel to streamlines; it acts 
on zero modes as the identity and commutes with /C, re- 
vealing a residual translation invariance for arbitrary 
that makes /C locally one-dimensional. Tj deforms states 
in directions perpendicular to streamlines: Tjf{J) = 
f{J + 77); it does not commute with IC in general, but 
does once K. is projected onto the eigenspace, reflect- 
ing an enhanced symmetry: the E — subspace may 
be decomposed into eigenstates of Tj and T^ simulta- 
neously. The deformations represent gapless modes, and 
are the counterpart of the area-preserving deformations 
that leave the mass limit of Hb invariant. The symme- 
tries of the continuum model are fully incorporated into 
the algebra of Poisson brackets, and a faithful discretiza- 
tion ought to preserve this algebra. 

V^ represents a singular perturbation to /C, and to 
proceed further /C must be regularized. On the lattice, 
the chain rule, which is needed to show that JCf{(j>) = 0, 
breaks down for local discretizations, so that a naive reg- 
ularization of JC destroys the conservation laws; at most 
a pair of eigenstates with zero energy survive. A simi- 



lar problem applies to a naive fourier discretization with 
momentum cutoff. Nevertheless, in taking the continuum 
limit, the energies of a subset of the modes approach zero; 
their eigenfunctions in the discrete theory are smoothly 
connected to those of the continuum theory as the UV 
cutoff is varied and their identification with zero modes 
of the continuum model remains. For this reason, the 
choice of regularization is not expected to affect physical 
quantities; however, if one needs to identify the subspace 
of zero modes explicitly, a symmetry preserving regular- 
ization is essential. 

The choice of regularization is therefore dictated by 
a need to preserve the algebraic properties of the con- 
tinuum model, even after truncating the system to a 
finite number of degrees of freedom. This require- 
ment is satisfied pOi by the Lie algebras su{N), where 



su{N) ^sdiff(T2) as N ^ oo. Setting p 
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in 0(x,p), one obtains the Heisenberg-Moyal bracket: 
[(l){x,-ih,dx),%l){x,-ihdx)]lh -^ {(l){x,p),'ip{x,p)}pB as 
h-^Q. Taking h — 2Tr/N for integer N yields, in fourier 
representation, the algebra su{N) with N'^ — 1 degrees 
of freedom. Note that this is the adjoint representation 
of su{N): ICijj = [(p, ijj], the wave-function is an operator, 
and the zero modes of /C are equivalent under unitary 
transformation to the A^ — 1 Casimir elements of su(N) 
together with the identity. 

In fourier space on the torus [0, l]x[0, 1] the Heisenberg- 
Moyal bracket takes the form: 

[0,V^]k = ^0k-i(27r/iV)sin((iV/27r)(k x l))Vi (3) 



where k, 1 = 27r(m, n)/N , < m,n < N ~ I. As shown 
in fig. la, for which V^(/) is taken as random and S- 
correlated on distances of order l/iV, this form correctly 
reproduces the overall spectrum obtained via lattice reg- 
ularization that is shown in fig. lb, except in the neigh- 
borhood of energy, where the N zero modes appear, 
separated from the rest of the spectrum by a pseudogap. 
Qualitatively, the gap can be thought of as originating in 
level repulsion; from eq. (2) it is clear that the low-lying 
states come from the neighborhood of saddle points of 0. 
The scaling of the size of the gap with N and its depen- 
dence on the correlations of cj) can be estimated by ob- 
serving that ii(j){k) ^ \k\^ then Vx(/)-V ~ Ar-('^+2); how- 
ever, 1/N of the spectral weight must be transferred to 
the zero modes, yielding a gap of order iV~(^+'^); £ = —2 
in fig. 1. 

My continuum arguments indicate that the eigenstates 
corresponding to Casimir elements of su{N) should be, 
in general, extended. I have confirmed this hypothesis 
by numerically evaluating their Chcrn number n^ . 

Chern number is computed by parameterizing the 
Hamiltonian by periodic variables 9^ , Oy and integrating 
the Berry flux over the Brillouin zone. The degeneracy 
of /C at energy leads to Berry matrices [O] rather than 
scalar Berry phases; to avoid this technical problem, I 
compute the Chern number for eigenstates of the oper- 
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FIG. 1. Density of states oi ilC in (a) su{N) discretization, 
and (b) N x N lattice discretization, averaged over many real- 
izations of (f). In (a) the A'^ states at _E = are suppressed for 
clarity. The inset expands (a) near _E = to exhibit the scal- 
ing of the pseudogap. Except near E — 0, almost all states 
are localized. 



ator Ti. = —V^ + ilC. There is no reason to expect that 
diffusion, or for that matter any other non-random, ho- 
mogeneous coupling, should localize otherwise delocal- 
ized states. It turns out that an extensive number of 
states with non-zero Chern number is found only within 
the 'diffusion band,' the energy scale set by — V^; these 
states span the entire diffusion band. 

Fig. 2 shows the fraction of eigenstates within the dif- 
fusion band having non-zero Chern number for a naive 
lattice discretization of ?i on an A^ x iV grid, as a func- 
tion of TV, averaged over many disorder realizations [p|. 
The components of A were chosen uniformly on an inte- 
val [—W, W] with (5-correlations in space, and the com- 
pressible part was subtracted. 9x and 6y represent phase 
twists applied at the boundary of the torus. The con- 
stant asymptote is consistent with a band of extended 
states in the thermodynamic limit. 

The inset shows Chern numbers for a typical disor- 
der realization, discretized with eq. (3) and N = 24. 
For this form of 7i, kinetic energy has been parame- 
terized as 2cos(27rm/iV + 9^,) + 2cos(27rn/Af + 9y) - 4, 
< 7TT.,n < A^ — 1. Nearly every eigenstate in the dif- 
fusion band has a non-zero Chern number; outside the 
diffusion band, only isolated energies, corresponding pre- 
sumably to the percolating streamlines, show non-zero 
Chern number. When i/C is replaced by V ^ all eigen- 
states have zero Chern number. These numerical results 
confirm extensive transfer matrix studies on the lattice 
@ for both RFI and PSI. 

The phenomena described in this paper originate in 
the infinite number of conserved fluxes that, for /C, are re- 
flected in a huge degeneracy at zero energy. The charges 
of the degeneracies produce non-zero Chern numbers 112] 
when z/C is perturbed by — V^, showing that the existence 
of extended states does not rely on the infinite degeneracy 
itself, but only on the fact that the perturbation doesn't 
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FIG. 2. Fraction / of states within diflusion band having 
non-zero Chern number, in lattice discretization of — V^ +iK.. 
The disorder amplitude W (see text) for circles 45.2; for 
squares, 22.6. Inset shows (integer) Chern numbers of eigen- 
states for a single realization of disorder and < -E < 16, in 
su[N) discretization. 



incompressible disorder, with its localizing effects most 
strongly felt near the band edges. For example, when A^ 
is negligible compared to — V^ + i/C then the entire diffu- 
sion band (or blob, for PSI) should be extended. When 
(j) ^ l/|k|, the critical energy ought to depend on the am- 
plitude of the vector potential randomess; as the relative 
strength of A^ increases, the extended regime should con- 
tract to the point of particle- hole symmetry . For these 
reasons, one might imagine that the divergent density of 
states recently found at the band center of a random-flux 
model 0] reflects the dirac sea I have identified above. 

I am grateful to P. Weichman, C. Mudry and B. Bialek 
for extensive comments on the manuscript, to D. Huse 
for observations about spatially anti-correlated flux dis- 
tributions, to M. Douglas for remarks on canonical trans- 
formations, and to R. Lifshits and M. Cross for support 
and encouragement. The numerics were carried out on 
the CACR parallel computer system at Caltech. 



lead to cancellation of the charges. 

More generally, one might expect that the infinite num- 
ber of fluxes, and the symmetry that generates them, 
would lead to an anomaly. This point of view applies nat- 
urally to vorticity w in a 2d incompressible fluid, where 
Jdaf{Lu) represents the constants of the motion corre- 
sponding to the /(0) discussed here. Vorticity pumped 
into a flow at large scales can be destroyed at small scales 
only by local action of V^, creating fluxes of conserved 
charge in momentum space. There is no reason to sin- 
gle out a particular charge, such as enstrophy, for special 
treatment, as in ref. p^ . 

I have exhibited an extensive fraction of delocalized 
states in the diffusion band when i/C 3> V^, that is, in the 
strong disorder limit, whereas arbitrarily weak potential 
disorder in d = 2 localizes all states of — V^ -I- V. Bar- 
ring some intermediate reentrant localized phase, which 
seems improbable, it is likely that the states will remain 
extended into the weak disorder regime iJC ^ V^. In 
RG language, one would naively expect fixed points at 
iK. = and DV^ = 0, both of which display an infinite 
number of extended eigenstates. 

The iK, <C V^ regime could be achieved if the random 
magnetic field B had strong anti-ferromagnetic correla- 
tions so that (j) r^ 1 and consequently A ^ fc. The center 
of the band sets a momentum scale kf, and the order 
of magnitude of the terms in RF may be estimated as 
V^ ~ fc?, i/C '^ kfk, and A^ ^ fc^, suggesting that it may 
be reasonable to incorporate the effect of the compress- 
ible disorder A^ as a perturbation to RFI. 

I conclude by addressing the effect of compressible dis- 
order on RFI. While eigenstates of V^ with opposite 
momentum are degenerate, as of course are modes of 
j/C, these degeneracies are broken for the extended eigen- 
states of — V^ -I- iJC; consequently, those eigenstates are 
chiral and small compressible disorder should not local- 
ize them. Rather, it can be expected to compete with 
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